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We consider the production of neutralino, chargino, and charged Higgs boson pairs in the MSSM 
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It is a widespread hope within the members of the high 
energy elementary particle physics community that the 
upcoming experiments at Tevatron [l| and at LHC |2| 
will finally reveal the existence of supersymmetric parti- 
cles, via direct production of sparticle-antisparticle pairs. 
In such an exciting hypothesis, the simplest available the- 
oretical supersymmetric model, the MSSM, will acquire 
the same role that belonged to the SM after the W, Z dis- 
coveries, and a long and crucial period of precision tests 
will start, aiming to confirm, or to disprove, the main 
theoretical assumptions that were used in the practical 
construction of the model. 

For the MSSM case, this program requires the con- 
struction of a linear e^e~ collider j^] that explores the 
cm. TeV energy region. This machine should be suffi- 
ciently accurate to provide the same kind of consistency 
tests of the model that were achieved in the hundred 
GeV region at LEP for the Standard Model, from de- 
tailed analyses of several independent one-loop virtual 
effects. 

If such an ambitious project is to be carried on, a 
clean investigation strategy will undoubtedly be welcome, 
given the fact that the involved theoretical model has a 
not simple structure in which several independent pa- 
rameters are implied. As a consequence of this feature, 
it might be appreciable to identify clean experimental 
measurements that, possibly, isolate in a relatively sim- 
ple way a reduced subset of parameters and of theoretical 
assumptions behind them. 

The aim of this short paper is to show that a ded- 
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icated combined analysis of the process of production 
of neutralino, chargino and charged Higgs boson pairs 
might lead to a relevant precision test of the theoretical 
assumptions that were used to fix, not only the basic su- 
persymmetric sector (i.e. the chiral superpotential and 
the gauge multiplet) , but also the gaugino component of 
the soft supersymmetry breaking of the model. This will 
be done in the paper by assuming a scenario of "mod- 
erately" light SUSY particles (i.e. the various masses 
are all below, roughly, 350 GcV) and a cm. energy y/s 
in the one TeV region. Under these two assumptions, a 
simple asyrnptotic expansion of so called "Sudakov type" 
is effective 4']. But the approach would not change for 
larger SUSY masses < 700 GeV since a related cm. en- 
ergy rescaling would still be possible, from our previous 
experience, up to ^/s values of about 2 TeV, where the 
resummation to higher orders is still not necessary |j| . 

After this brief description of the strategy of our ap- 
proach we are now ready to illustrate the practical de- 
tails of our work. With this aim, we must recall, for the 
reader's convenience, that the relevant details of an anal- 
ogous procedure, limited to the combination of chargino 
and charged Higgs pairs production, have already ap- 
peared in a previous paper |5|- It was shown there that 
from the combined analyses of the slopes of the total cross 
sections in a cm. energy range of about 1 TeV, assuming 
a "light" SUSY scenario where all the relevant sparticle 
masses lie below ~ 350-400 GeV , a strong constraint on 
tan /3 (mostly produced by the Higgs cross section) and a 
strip in the (M2 , /.t) plane (only produced by the chargino 
cross sections) were derivable. Here M2 is the soft SUSY 
breaking wino mass, ^ is the Higgsino mass parameter, 
and tan/3 — ^2/^1 is the ratio of the two Higgs doublets 



vevs. The novel process that we consider in this paper 
is neutrahno pair production, and we shall show that its 
addition to the other processes will lead to highly im- 
proved constraints on the previous parameters, with the 
additional presence of Mi, the soft SUSY breaking bino 
mass. With this aim, we now briefly list the relevant 
theoretical formulae that we need for the analysis. 

At the Born level, the neutralino pair production pro- 
cess e+e^ -^ Xi'Xj is described by three s, t, u channel 
components of the scattering amplitude, using the fol- 
lowing notations (the indices a, b label the electron and 
the z-th neutralino chiralities) : 
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We assume a previous production of the charginos, of the 
charged Higgs bosons and of at least two neutralinos with 
mass Mo , Mo . Calling AfgusY the heaviest of the real 
and virtual SUSY particles that appear in the processes, 
and assuming a "reasonable" limit Msusy — 350 — 400 
GeV, we shall choose the value -/s ~ 1 TeV value to pro- 
ceed with our approach. This is due to the fact that, in 
a previous paper only concerned with the charged Higgs 
pair production [2|, we proved that in such a configura- 
tion an asymptotic energy logarithmic expansion of so 
called "Sudakov type" was reliable, with the only ad- 
dition of a constant term to the leading quadratic and 
next-to-leading linear logarithmic terms. This conclusion 
allows to propose a determination of the SUSY parame- 
ters entering the Sudakov logarithms, based on measure- 
ments of the slope of the total cross section, in which the 
(complicated) constant terms cancel. 



The approach that we shall follow in this paper, that 
was already used in the combined chargino-charged Higgs 
analysis [5J, will assume that a similar expansion is valid, 
with only logarithmic and constant terms. To prove 
this statement would require a detailed comparison of 
the complete existing calculation jg with the assumed 
asymptotic expansion. This will be the goal of a forth- 
coming rigorous analysis. For the moment, we shall as- 
sume its validity as a working Ansatz, and we shall show 
the main relevant consequences that it will be able to 
produce. 
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After this preliminary discussion, we are now ready to 
(6) write the relevant asymptotic expansion of the scattering 
amplitude at one loop, in which we shall only retain the 
leading quadratic and next-to-lcading linear logarithmic 
terms. For this purpose, we shall use the following for- 
(' ) mal decomposition the separate F — S, T or U subam- 
plitudes (the Born values are those listed in Eqs. H2I8(I ^: 
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The quantities Z^- are the elements of the 4x4 mixing 
matrix, defined in a conventional way |j] . For simplicity, 
we have neglected in the considered asymptotic region 
the selectron mass in Eq.(l). 

For the purposes of a precision test, it becomes manda- 
tory to compute the following perturbative one-loop ex- 
pansion of the scattering amplitude. This requires the 
calculation of a large number of Feynman diagrams. The 
resulting expressions are valid for any value of the cm. 
energy y/s = \/{Pe- +Pe+1^- The full one-loop result has 
been calculated recently [y|. In this work we propose an 
alternative approach, valid for values of y/s "much" larger 
than all the SUSY sparticle masses involved in the pro- 
cess. In fact, our strategy should be now made very clear. 
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In the asymptotic expansion, two different kinds of log- 
arithmic terms appear. The first ones are the standard 
linear ones of RG origin. They are well known and can 
be derived in a straightforward way replacing in the Born 
quantities the various bare couplings with running ones. 
For sake of completeness, we write the related formulae, 
noticing that they are only requested for the s-channel 



(purely Higgsino) component, as: 
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where (3o = —1/4 and /3q = —11/4. The second type 
of logarithms which arise asymptotically is that of the 
"genuine weak" Sudakov terms. These are usually classi- 
fied [8j as logarithms of gauge non universal, gauge uni- 
versal and Yukawa origin. The gauge non universal, scat- 
tering angle dependent ones, stem from box diagrams and 
t,u channel vertices and have the following expressions: 
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The gauge universal ones are due both to the initial and 
to the final vertices. The initial contribution is fixed by 
the coefficients 
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The contribution from final neutralino legs is different for 
the s or u, t channels. For the s-channel we have higgsino 
components 
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The logarithms of Yukawa origin are only due to final 
s-channel vertices with virtual heavy (b,t) quark-squark 
pairs. These are the only logarithmic terms that contain 
(also) the vevs ratio tan/3 and their expression is 
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To conclude, the final contribution to the u and t channel 
amplitudes is due to gaugino components for which there 
is no Yukawa contribution: 
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Eas. (|10l25|l represent the complete logarithmic contri- 
butions at one loop to the considered process, and are the 
main original result of this paper. We expect from our 
previous discussion their reliability for what concerns the 
calculation of the slope of the total cross section in the 
1 TeV range. To perform the latter, one must first com- 
pute the expression of the differential cross section that is 
readily obtained from the above amplitude. After angu- 
lar integration, one obtains the approximate (to next-to- 
leading logarithmic order) asymptotic expression of the 
total cross section. This will be sufficient to determine 
the effective asymptotic expression of the slope of the 
cross section, where by definition possible extra "next- 
to-next-to-leading" (i.e. constant) terms disappear. This 



TABLE I: Input parameters and masses of charginos, and 
lightest neutralinos for the three input sets Si, S2, and S3. 
The masses are expressed in GeV. 



RR2: M| = 78 GeV, M, = 150 GeV, |l = 263 GeV, Ian |3 = 30 
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expression will only depend on tan [3 and on the mixing 
parameters Zf, . The latter ones, in turn, can be ex- 
pressed as functions of the supersymmetric parameters 
Ml, M2, fx and tan/3. Thus, the final form of the slope 
will depend on these four parameters. 

In order to achieve the maximal theoretical informa- 
tion using our approach, we have combined the theoret- 
ical expressions of the slope of the neutralino pair to- 
tal cross section, given in this paper, with those of the 
chargino and charged Higgs pairs, given in pj. In prac- 
tice, we have limited our analysis to the production of 
the two Hght observable final states {xiX2' X2X2) fo'^ '^su- 
tralinos, combined with the production of all the three 
chargino pairs and of the charged Higgs pair. To work 
in a consistent way, we have chosen as relevant examples 
for the masses of the produced pairs three sets of values 
denoted by S'l, 5*2, and 5*3. The first, Si, is the Tesla 
benchmark point RR2|9j with the two lightest neutrali- 
nos being respectively 95 % bino and 82 % wino. The 
set ^2 is a mixed scenario with neutralinos having non 
negligible gaugino and Higgsino components; Xi is 86 % 
bino and 13 % higgsino, X2 is H % bino, 48 % wino and 
41 % higgsino. Finally, 6*3 is a purely Higgsino one with 
the two lightest neutralinos being 92 % and 98 % hig- 
gsino like. The values of the input parameters as well as 
the masses of the two charginos and of the two lightest 
neutralinos are summarized in Tab. (1). We have per- 
formed a standard y^ analysis assuming in the various 
scenarios 10-12 experimental points at energies ranging 
from 700-850 GeV (depending on the scenario) up to 1200 
GeV and assuming an experimental accuracy of 1% for 
the cross sections of gauginos and 2% for that of Higgs 
bosons. For each scenario we have checked the condition 
y/s ^ AfgusY- The results of our analysis are shown in 
Figs. (1-3). From inspection of the figures one can draw 
the following conclusions. In Fig.(l) we consider the 
S'l benchmark point. One notices that a combined anal- 
ysis is able to generate closed contours in the planes of 
all possible six couples of parameters. This possibility 
is obviously existing for the first 3 cases i.e. (Mi,M2), 
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FIG. 1: Si benchmark point, la error bounds on the MSSM 
parameters M2, /i, tan/3. In this and in the following figures 
the crosses denote the values of the parameters in the specific 
benchmark. 



(Mi,^), (Ml, tan /3) under the condition that the neu- 
tralino information is added to the remaining ones, that 
do not depend on Mi . Less obvious and more illustrative 
is the closure of the undetermined strip in the (Af2,M) 
plane (drawn in the absence of neutralino data) when 
neutralino data are added to the remaining ones. This 
addition does not practically improve, on the contrary, 
the bounds on (M2,tan/3) and (/i,tan/3) obtained from 
chargino-Higgs data. One sees typical errors, under the 
assumed experimental conditions, of about 10 GeV for 
Ml, M2, II and of a relative -- 30 - 40 % for tan/3. In 



Mixed: M, = 100 GeV, M, = 200 GeV, n = 200 GeV, Ian |3 = 30 
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FIG. 2: S2 benchmark point, la error bounds on the MSSM 
parameters M2, pi, tan/3. 

Fig. (2) we consider the ^2 scenario. One finds, approxi- 
mately, the same features for the set of projections. The 



size of errors is now of approximately ^ 10 GeV for Mi, 



Higgsino scenario: M = 200 GeV, M, = 400 GeV, n = 100 GeV, Ian p = 30 
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FIG. 3: S3 scenario, la error bounds on the MSSM parame- 
ters M2, /i, tan/3. 

S3, a typical Higgsino scenario. Here the new (negative) 
feature concerns the lack of boundary for Afi, as seen 
from the first three plots. This can be understood since, 
at the considered Mi, M2 and /i values, the mixing pa- 
rameters Zj^ do not vary appreciably when Afi increases. 
On the contrary, the three plots in the bottom line retain 
the typical features of the previous cases, with particular 
relevance of the neutralino role for the determination of 
lower limits for fi and M 2 . The errors are now , typically, 
of about 20 GeV for M2, /i and of a relative ~ 10% for 
tan/3. 

The previous results should illustrate the aimed possi- 
ble outcomes of the testing strategy that we are proposing 
in this preliminary qualitative paper. Our assumed phys- 
ical inputs will necessarily be the masses of the (suppos- 
edly produced) (light) neutralinos, charginos and charged 
Higgs. Their values will depend on the four parameters 



Ml, M2, n, tan/3 but, in general there will not be a 1-1 
correspondence, and different sets of parameters might 
reproduce "essentially" the same masses [llj. Within 
the errors that we have illustrated in our examples, these 
different sets can be discriminated via the x^ analysis 
that we have proposed, and the "correct" set is selected 
by the true experimental data, if the latter are actually 
described by the model. If none of the candidate sets 
were compatible with the analysis, an indication would 
arise that some of the details of the model might need a 
suitable modification. 
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